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In terms of a new QCD factorization formula for ^/ip production, we calculate the ^/ip 
suppression in nuclear collisions by including the multiple scattering between the pre- 
J/V" partonic states and the nuclear medium. We find agreement with all data on ^/ip 
suppression in hadron-nucleus and nucleus-nucleus collisions, except a couple of points 
(the "second drop" ) at the highest Et bins of the new NA50 data. 

1. J/^ Production 

^/ip suppression in relativistic heavy ion collisions was suggested as a potential signal of 
the quark-gluon plasma [|I|. In order to study the i /ip suppression quantitatively, we need 
to understand how the i /ip mesons are produced, and the reliability of QCD calculations 
for i/ip production. 

In order to produce a. ^/ip meson in hadronic collisions, the energy exchanges in the 
collisions have to be larger than the ^/ip mass (Mj/^) or the invariant mass of the produced 
quark pair {rricc > 2Mc ~ 3 GeV). The cc pairs should be produced at a distance scale 
th < l/2Mc ~ 1/15 fm. Since this is much smaller than the size oi a, ^/ip and the energy 
exchanges are much larger than the nonperturbative momentum scales in the ^/ip wave 
functions, the J/?/' is unlikely to be formed at the collision point. Instead, the ^/tp meson 
should be formed after some resonance interactions (or coherent soft gluon interactions 
between the charm and anticharm quarks). Thus, the transformation from a pre- J/'?/' 
cc pair to a physical J/^ occurs over several Fermi ||^. During this transformation, the 
produced cc pairs can also radiate gluons and have interactions with spectators. Due to 
gluon radiation and the spectator interactions, it is possible for the cc states, even with 
different quantum numbers, to become a J/ip meson. 

However, because of the large energy exchange, the spectators' roles in the hard colli- 
sions are suppressed by an extra factor piyL/[{2McY + q^] with the medium density 
and medium length L. If the transverse momentum qt is large enough, the spectator 
interactions are strongly suppressed, and the cross sections for hadronic J/tp production 
can be factorized 0, 



daAB^J/ip 

dyd'^qT 



/ dx(pa/A{x) / dx'(pb/B{x') 



a,b 




(1) 



where J2ab ^^un over all parton flavors and the (p{x) are the normal parton distributions. 
The (Tafe^j/^(a;, x') represents the cross section for two partons to produce a J/ip. Without 
the nuclear medium, L is of the order of the hadron radius (i? ~ 1 fm). 
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Because the production of the pie-J/tp cc states and the transformation to the J ftp take 
place at very different distance scales, quantum interference between these two stages is 
suppressed by the ratio of these two distance scales. Therefore, the partonic cross sections 
for the J/ip production in Eq. (|^) can be further factorized into two stages: (1) short- 
distance production of the pre-J/ip partonic states, and (2) long-distance transition from 
the partonic states to the physical J/ip mesons 0, 

d^ab^j/^ sr^ f dz f 2 1 da ah 



dyd^qT 'riJ J \dm'l^dyd?p, 



^ / — / dml, Dy,^^„^{z- ml,) (2) 



cct , 



where X][cc] sums over all cc states, pcc = Pij^jz and m^g -C |pccP- The D\c^ represent 
transition probabilities from pre-J/?/^ cc states to J/V' mesons, and are defined by matrix 
elements of non-local operators 0. The differences between our approach and other 
models of J/V' production can be cast as differences in -D[cc]- 

In QCD perturbation theory, the Dyg^ include all soft resonance interactions and gluon 
radiations. Similar to the parton-to-hadron fragmentation functions, the Z)[cc] depend on 
nonperturbative physics. Thus, the predictive power of the factorized formula relies on 
the universality of these transition distributions 

The non-relativistic QCD (NRQCD) model 0] and the color evaporation (CE) model 
P correspond to two diverse assumptions for the D^cc\- The NRQCD model assumes that 
the D[cc] are steeply falling functions of the invariant masses of the states, and sensitive 
to the pairs' quantum numbers (such as color and spin). We can represent the leading 
contributions of this model by expanding the partonic cross section aah^\c^ in Eq- (@) 
at m^g = 4M^ and z = \. In the NRQCD model, the produced cc pairs with large m^g 
have very small probabilities to become J/?/' mesons. On the other hand, the CE model 
assumes that the D[cc] are independent of the pairs' invariant masses (provided they are 
below the open charm threshold), and are not sensitive to the pairs' quantum numbers. 
We can represent the CE model by taking Dycc\ = fj/tp^i^ — z)9{4:M^ — m^g) in Eq. (||) 
with a fitting constant /j/^ independent of the cc pairs' quantum numbers. 

Since the NRQCD and the CE models are both consistent with the existing data on 
inclusive J/ip production, the new QCD factorized formula should also be consistent with 
the data [^]. However, both models have difficulties explaining the recent CDF data from 
Fermilab on J ftp polarization at large transverse momentum |P. 

At large transverse momentum, J/ijj are produced mainly from the fragmentation con- 
tribution, as shown in Fig. 1, where a virtual gluon is produced and decays into a cc pair 
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which transmutes into a physical i/ip. The polarization of the produced i/ip depends 
on the polarization of the cc pair and the details of the transition distributions -D[cc]. In 
order to use the polarization measurements to discern between the production models, we 
need independent tests of the reliability of the QCD formalism for calculating the cc po- 
larization at large transverse momentum. As illustrated in Fig. 1, the low mass Drell-Yan 
lepton-pair's angular distributions at large transverse momentum have a lot in common 
with the cc polarization 0. If the QCD calculations are consistent with the virtual pho- 
ton polarization in Drell-Yan massive lepton-pair production, we can conclude that the 
QCD formalism for calculating the polarization of the cc at large transverse momentum is 
reliable, and therefore, measurements oi^/ip polarization provide decisive tests for models 
of J/V' production. 

2. J/^ Suppression 

In addition to polarization measurements, the nuclear medium dependence of J/?/) pro- 
duction is sensitive to the ^/ip formation process [§]. If the ^/ip mesons are formed 
immediately after the production of the cc pair, the medium dependence is dominated by 
the hadronic scattering between a color-singlet ^/ip and the nuclear matter. On the other 
hand, if the ^/ip mesons are formed much later, the observed medium dependence is a 
consequence of the dynamical interactions between the pre-J/V' cc states and the nuclear 
medium. Therefore, the nuclear dependence can help select the correct model for J/^ 
production. 

In proton-nucleus or nucleus-nucleus collisions, the produced cc pairs are likely to inter- 
act with the nuclear medium before they exit. If we assume each interaction between the 
cc pair and the nuclear medium is about the same and can be treated independently, we 
naturally derive the Glauber formula: aAB = PN(^ahsL{A, B)] with effective absorp- 

tion cross section cXabs to change a cc pair to a pair of open charms, and effective length 
L{A, B) of the medium in A+B collisions. The characteristic feature of the Glauber for- 
mula is a straight line on a semi-log plot of aAB vs. L{A, B), which is not consistent with 
the strong J/ip suppression observed in Pb-Pb collisions ||^. 

Our model of J/ip suppression is based on the following points: (1) it is the cc pair, not 
the J/ip meson, that interacts through most of the nuclear medium; (2) multiple interac- 
tions of the c and/or c with soft gluons in the medium increase the relative momentum 
between the c and c; and (3) as the relative momentum increases, the phase space to 
form a J/tp meson decreases |p. Since multiple soft-scatterings lead to a larger relative 
momentum between the c and c, the effective (Tabs of the pair increases as the cc passes 
through the medium leading to a stronger suppression than given by the simple Glauber 
formula which applies to single particle propagation p . 

Integrating over J/ip^s momentum in Eq. (|T]), we obtain a factorized expression for the 
J/ip total cross section 

crAB-~^.yi, ~ H / dx(j)a/Aix) J dx' (pb/ b{,x') j dml^^^^^^^^F[cc]^,y^{ml^) (3) 

a,b,[cc] cc 

where F[cc]^3/,p represents an inclusive transition probability for a cc pair of invariant 
mass to become a. J/ip meson. Without the large transverse momentum, the cor- 
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rections to Eq. (^) from the spectator interactions, which are suppressed by a factor 
~ pi\iL{A, B)/AM^, become more important because of the large medium length L{A, B). 
Just like the random walk of two particles, the soft spectator interactions effectively in- 
crease the invariant masses of the cc pairs, and therefore, reduce the phase space to form 
mesons. We calculate the ^/ip suppression in the nuclear medium by including the 
multiple spectator interactions to the pre-J/V' cc pairs, and plot our results in Figs. 2 and 
3. The data are from Ref. P|. The dashed lines are predictions of the Glauber formula, 
and the solid lines are the results of our calculations. With only one parameter e^, which 
is defined as average gain of invariant mass per unit medium length and is fixed by the 
total suppression in Fig. 2, our calculations of the i /ip suppression are consistent with all 
existing data, except a couple of points at the highest Et bins of the New NA50 data [§]. 

We conclude that the NRQCD model and the CE model are closely connected to the new 
QCD factorization formula for ^/ip production. The measurements oi ^/tp polarization 
could provide a decisive test of the production mechanism. Our model of J/?/) suppression 
is consistent with nearly all existing data on ^/ip total cross section. 

This work was supported in part by the U.S. Department of Energy under Grant No. 
DE-FG02-87ER40731. 
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